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We propose a generic test for models in which gravity is modified to do away with dark matter. 
These models tend to have gravitons couple to a different metric than ordinary matter. A strong 
test of such models comes from comparing the arrival time of the gravitational wave pulse from a 
cosmological event such as a supernova with the arrival times of the associated pulses of neutrinos 
and photons. For SN 1987a we show that the gravity wave would have arrived 5.3 days after the 
neutrino pulse. 

PACS numbers: 04.80.Cc,95.85.Sz,95.85.Pw,95.35.+d 



1. Introduction 

As early as 1933 Zwicky was able to infer the need for 
dark matter by applying the virial theorem to observa- 
tions of the Coma Cluster In the 1970's Rubin, Ford 
and Thonnard accumulated independent evidence from 
the rotation curves of spiral galaxies 0, 01, 0] • And the 
1990's saw weak lensingused to probe dark matter in 
galactic clusters @, H, [ill, H, [Io| . 

It will be seen that the evidence for dark matter has 
so far been entirely restricted to gravity. That is, one 
infers the Einstein tensor from measurements of cosmic 
motions, or from lensing, and then compares the result 
with the observed stress-energy in stars and gas, 
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There is indeed a disagreement between the two sides of 
the equation, but it is not clear that this signifies the 
need for dark matter on the right of U) rather than a 
modification of gravity on the left. Although there are 
plausible dark matter candidates, none of them has so 
far been detected in a terrestrial laboratory 11, 12, 

Certain regularities in cosmic structures suggest mod- 
ified gravity. One of these is the Tully-Fisher relation, 
which states that the luminosity of a spiral galaxy is pro- 
portional to the fourth power of the peak velocity in its 
rotation curve If luminous matter is insignificant 



compared to dark matter, why should such a relation ex- 
ist? Another regularity is Milgrom's Law, which states 
that the need for dark matter occurs at gravitational 
accelerations of ao ~ 10~ 10 m/s 2 (H)]. This has been 
observed in cosmic structures which vary in size by six 
orders of magnitude (lfjj . 

A modification of Newtonian gravity which explains 
these regularities was proposed by Milgrom in 1983 [l7j |. 
His model, Modified Newtonian Dynamics (MOND), was 
soon given a Lagrangian formulation in which conserva- 
tion of energy, 3-momentum and angular momentum are 



manifest [l8j. However, there was for years no success- 
ful relativistic generalization which could be employed 
to study cosmological evolution. Even in the context of 
static, spherically symmetric geometries, 
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the early formulation of MOND fixed only B(r), not 
A(r). It was therefore incapable of making definitive pre- 
dictions about gravitational lensing. 

A relativistic extension of MOND has recently been 
proposed by Bekenstein [llj. This model is known as 
TeVeS for "Tensor- Vector-Scalar." In addition to repro- 
ducing the MOND force law at low accelerations, TeVeS 
has acceptable post Newtonian parameters, and it gives 
a plausible amount of gravitational lensing [lj^. When 
Te VeS is used in place of general relativity + dark mat- 
ter to study cosmological evolution, the results are in 
better agreement with data than many thought possible 

H ED, 0,(23. 

What concerns us here is the curious property of TeVeS 
that small amplitude gravity waves are governed by gen- 
eral relativity through the metric g^ v , whereas matter 
couples to a "disformally transformed" metric which in- 
volves the vector and scalar fields, 
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The Scalar- Vector- Tensor gravity (SVTG) theory pro- 
posed by Moffat also has different metrics for matter and 
small amplitude gravity waves 0, . The appearance 
of this feature in two very different models is the result of 
trying to reconcile solar system tests with modified grav- 
ity at ultra-low accelerations. Solar system tests strongly 
predispose the Lagrangian to possess an Einstein- Hilbert 
term [26[. On the other hand, failed attempts to gener- 
alize MOND 27J have led to a theorem that one cannot 
get sufficient weak lensing from a stable, covariant and 
purely metric theory which reproduces the Tully-Fisher 
relation without dark matter [28JI. Hence the MOND 
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force must be carried by some other field, and it is a 
combination of this other field and the metric which de- 
termines the geodesies for ordinary matter. However, the 
dynamics of small amplitude gravity waves are still set by 
the linearized Einstein equation. This simple observation 
makes for a sensitive and generic test. 

We define a Dark Matter Emulator as any modified 
gravity theory for which: 

1. Ordinary matter couples to the metric g^ that 
would be produced by general relativity + dark 
matter; and 

2. Small amplitude gravity waves couple to the met- 
ric produced by general relativity without dark 
matter. 

Now consider a cosmic event such as a supernova which 
emits simultaneous pulses of gravity waves and either 
neutrinos or photons. If physics is described by a dark 
matter emulator then the pulse of gravity waves will 
reach us on a lightlike geodesic of g^, whereas neutri- 
nos and photons travel along a lightlike geodesic of g pv . 
If significant propagation occurs over regions that would 
be dark matter dominated in general relativity then there 
will be a measurable lag between arrival times. 

The magnitude of the expected time lag is so large 
that great precision is not needed for either g^ or 'g^y. 
We demonstrate this in section 2 by computing the time 
lag for SN 1987a using simple models for the distribu- 
tion of luminous and dark matter. Section 3 discusses 
the prospects for observation. Our conclusions comprise 
section 4. 



2. Time Lag 

We model the luminous matter of our galaxy as a 
spherical bulge of M ~ 10 41 kg. For this source general 
relativity gives a metric of the form ([2|), with the func- 
tions B and A expressed in terms of twice the Newtonian 
potential e = (2GA/)/(c 2 r), 
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B(r) =l-e , A{r) = — 
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At the Sun's radius of rs — 8.0 kpc we have e ~ 6 x 10~ 7 , 
so it is valid to neglect terms of order e 2 . 

To determine ~g^ v we assume the galaxy is surrounded 
by an isothermal halo which makes the asymptotic rota- 
tion speed i>* = (aoGM)^. This gives rise to a second 
small parameter, e* = 2v 2 /c 2 ~ 6 x 1CP 7 . For this source 
general relativity gives a metric of the form ((2]) with, 
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We have chosen the integration constant in the order e* 
contribution to B{r) so that it vanishes at r — rs- In 
section 4 we discuss the effect of other choices. 



It is most efficient to work this problem in Cartesian 
coordinates for which the nonzero components of the two 
metrics are, 
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Here the radial unit vector is r l = x l /r, the radius is 
r = 1 1 a; 1 1 and we have neglected terms which are higher 
than linear in either e ore,. 

In each metric we must construct the lightlike geodesic 
from a point xl in the Large Magellanic Cloud to the 
Sun's position x*s- For g^ this means solving the second 
order equation, 
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The general solution to jS} depends upon eight integra- 
tion constants which can be taken as the initial position 
X M (0) and velocity X M (0). Because g^ v is static we may 
as well start the geodesic at t = 0, 
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(9) 



We fix the spatial components of the initial velocities by 
requiring the geodesic to reach the Sun at r = 1, 



X l (l)=^s 
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The initial temporal velocity is fixed by requiring the 
geodesic to be lightlike, 



g^(0,x L )r(0)r(0) = 
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The geodesic x M (t) is determined the same way, but with 
the substitutions of for T^ pa . in |J5]) and of g^ for 
g^v in (fTTjl . The time lag we seek is, 
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Because e and e* are so small we can use perturbation 
theory. The th order (flat space) solution is the same 
for both geodesies, 



X°(t) = Axt 



(13) 



where Ax* = x % s — x % L . Because both g^ and g^ have 
the same order e terms, the order e corrections to 
cancel those to x p in (|12[) and we need only consider the 
corrections of order e*. The nonzero components of the 
order e* connection are, 
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The 1 st order corrected geodesic follows from integrating 
the geodesic equation, 

Ax fl (r)sf(r)-X"(r), (15) 
= A*"(0)- fair' fair" Ar» pa (x(r"))ifx ■ (16) 
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The spatial boundary conditions ([TO]) imply, 



Ax'(0) = /dr(l-r)Ar' x(r) W$ 



(17) 



The lightlike condition (fTTjl fixes the final integration con- 
stant, 



AX°(0) = A -'' A ^ :lli ' A 
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The parameter integrations are straightforward when ex- 
pressed in terms of the dimensionless constants, 



a = 



xi ■ Ax 



and 



Ax 2 ' Ax 2 

Our result for the first order time lag is, 

e»Aa; 
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Within the context of our simple model for the mat- 
ter distributions, equation (|20p is valid for any xl- The 
specific values appropriate to SN 1987a are [291 ], 

r s ~ 8.0 kpc , r L ~ 50.9 kpc , Ax ~ 51.4 kpc , 
a ~ -.9775 and (3 ~ .9793 . (21) 

Substituting these values in (|20|) gives, 
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36.7 days 



.144 



-5.3 days 
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Therefore simultaneous pulses of neutrinos and/or pho- 
tons would have arrived 5.3 days before the gravity waves 
if physics were described by a dark matter emulator. 

The sign of ([2H)l is mostly set by the difference of the 
two gravitational potentials, B(r) — B(r) ~ e* ln(r/rs)- 
The lag is negative because ordinary matter is at higher 
potential than gravity waves. The square bracketed term 
in (J20J) is typically of order one. (It is -.144 for SN 
1987a.) Hence the rough magnitude of the time lag is 
e* ~ 6 x 10 -7 times the fiat space propagation time. 
This works out to be 36.7 days for SN 1987a. One could 
therefore anticipate a lag in the range of days without 
the detailed computation. 



matter. If physics is described by general relativity + 
dark matter then the arrival times will be simultaneous. 
If physics is instead described by a dark matter emulator 
then the photon and neutrino pulses will generally arrive 
some days before the gravity wave signal. 

One only expects a few supernovae per century in our 
galaxy. Of course the light signal from a supernova can 
be observed even from distant galaxies, but this is much 
further than we have any hope of detecting neutrinos or 
gravity waves. The optical signal from a supernova will 
typically lag the neutrino pulse by several hours — as 
was the case for SN 1987a [29( — because photons must 
traverse the optically dense stellar environment. How- 
ever, the magnitude of the expected delay for the gravity 
wave signal is so great that this should not matter except 
for nearby supernovae. 

The amount of gravitational radiation from a super- 
nova depends on the oblateness of the progenitor. If the 
oblateness of SN 1987a was in relation to that of the Sun 
then current gravitational wave detectors would prob- 
ably not have seen anything (30j . However, advanced 
LIGO would detect such a supernova out to 0.8 Mpc 
[30j ] . This includes the Andromeda galaxy, which dou- 
bles the expected rate and also ensures that the signal 
passes through dark matter dominated regions. 

Neutrinos from SN 1987a were observed by the Kami- 
okandc-II [3l|, [32j and Irvine-Michigan-Brookhaven 33j, 
HH detectors. Of course the state of the art has im- 
proved since then, although effective coverage is still lim- 
ited to our own galaxy and its satellites 1351. 1361 . There is 
now a Supernova Early Warning System [37[ comprised 
of Super-Kamiokande [38|], the Large Volume Detector 
[3^ |. the Sudbury Neutrino Observatory — which is 
no longer operating but may be succeeded by SNO+ 41 1 
— and the IceCube Detector 1 121. 



4. Discussion 

Our largest uncertainty is modeling the radius at which 
the density of dark matter reaches its asymptotic form of 
p(r) ~ v 2 /{A-KGr 2 ). We chose that radius to be rs in ([5|) 
because this is near where the Newtonian acceleration 
reaches the value of ao ~ 10~ 10 m/s 2 at which Milgrom's 
Law implies the need for dark matter 15[ . Had we chosen 
some other radius, rp, the square bracketed term in (|20|) 
would change to 



3. Observational Prospects 

The test we have proposed requires detecting gravity 
waves and either neutrinos or photons that were emitted 
simultaneously from some cosmic event separated from us 
by an extensive region in which the gravitational influ- 
ence of dark matter predominates over that of luminous 
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One would expect this to make the time lag more negative 
because it is difficult to imagine having dark matter "turn 
on" much beyond rs- 

If neutrinos have mass then a neutrino of energy E 
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would propagate with velocity, 



E 2 



(24) 



Taking the average detected supernova neutrino energy 
to be E ~ 10 MeV [Zj|, and the (electron) neutrino mass 
to be its current upper limit of m ~ 2 eV [44j , the velocity 



would be changed by Av ~ — 2 x 10 



Over a distance 



of 50 kpc this would delay the neutrino by .1 s, which is 
negligible on the scale of days expected for the primary 
effect. 

The generality of our test derives from avoiding the 
details of specific models and concentrating instead upon 
what they do, which is to emulate dark matter. Rota- 
tion curves and lensing are too well measured for any 
viable model not to somehow reproduce the metric 
which couples to ordinary matter. We have explained the 
theoretical reasons why gravity waves are predisposed to 
couple to the original metric g^ v [2ri I2II . Both of the 
existing models do have this property [lij 0] , but as we 
cannot prove it to be unavoidable, it is safest to comment 
that our results are limited to dark matter emulators. 
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